Local Enhancement of Antiferromagnetic Correlations by Nonmagnetic Impurities 
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The local enhancement of antiferromagnetic correlations near vacancies observed in a variety of 
spin systems is analyzed in a single framework. Variational calculations suggest that the resonating- 
valence-bond character of the spin correlations at short distances is responsible for the enhancement. 
Numerical results for uniform spin chains, with and without frustration, dimerized chains, ladders, 
and two dimensional clusters are in agreement with our conjecture. This short distance phenomenon 
occurs independently of the long distance behavior of the spin correlations in the undoped system. 
Experimental predictions for a variety of compounds are briefly discussed. 
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Studies of ladder compounds continue producing fas- 
cinating results. In addition to the discovery of a spin 
gap in undoped even- leg ladders jl[], superconductivity 
at high pressure in Sro.4Ca13.6Cu24O41.s4, with 2-leg lad- 
ders and chains in its structure, has been recently re- 
ported ||. Both properties, predicted by theoretical ar- 
guments, [|| indicate a close interplay between the spin 
and charge degrees of freedom leading to a rich phase 
diagram. More recently, the doping of ladders with non- 
magnetic impurities (replacing spin 1/2 Cu 2+ by spin 
Zn 2+ ) has revealed another surprising property: the 
spin gap is rapidly suppressed as the Zn concentration 
increases, and an antiferromagnetic (AF) phase is sta- 
bilized 0]. A similar behavior has also been observed in 
spin-Peierls chains , which have a spin gap produced by 
dimerization. The phenomenon is interesting since a spin 
ordered state is generated by the random replacement of 
spins by vacancies, an apparently disordering procedure. 
These results have been recently addressed with one di- 
mensional (ID) spin models using field theory and 
numerical techniques. Computational studies found that 
the AF correlations near a vacancy in dimerized chains 
and 2-leg ladders are enhanced with respect to 
the undoped case. It was conjectured that this local en- 
hancement may trigger the 3D AF order in Zn-doped 
dimerized chains and ladders. In-gap weakly interacting 
5=1/2 localized states were found near Zn How- 
ever, the microscopic origin of the local AF enhancement 
near a vacancy is still not intuitively understood. 

Independently of these recent developments, related 
phenomena have been discussed in a variety of contexts: 
1. A staggered moment appears near a vacancy for ID 
5 = 1 Heisenberg systems [[To); 2. The undimerized ID 
5=1/2 Heisenberg model has an enhanced spin struc- 
ture factor S(tt) near vacancies according to boundary 
conformal field theory and Monte Carlo (MC) simula- 
tions |llj; 3. Near a vacancy injected into a 2D Neel 
ordered state, the staggered moment increases with re- 



spect to the undoped system [ |12| . 

In this paper it is proposed that all these examples 
of locally enhanced antiferromagnetism near a vacancy, 
which have been studied independently in the literature, 
may have a simple common explanation. The unify- 
ing picture relies on the resonating- valence-bond (RVB) 
character of the spin correlations at short distances for 
a variety of Heisenberg spin systems where the nearest- 
neighbors (NN) interaction, regulated by a coupling J\, 
is the largest. Independently of the long distance proper- 
ties of the model, the small distance behavior, at least for 
small spin 5, is dominated by the formation of short spin 
singlets resonating in all their possible arrangements [jl3| . 
Actually, RVB variational states have large overlaps with 
exact ground states of finite clusters, even for 2D Neel or- 
dered systems 0] . Corrections involve long singlets and 
triplets with a weight decaying as the spin-spin distance 
grows. 

In the case of a ID 5 = 1/2 chain, Fig.la-b illustrates 
two arrangements of spin singlets that carry important 
weight in the undoped ground state. Consider now a va- 
cancy introduced at site 0, effectively cutting the chain 
if a NN spin Hamiltonian is used. In this situation con- 
figuration Fig. lb is replaced by Fig.lc that contains a 
"free" spin. The NN exchange strongly favors the cou- 
pling of this free spin with the spin at site 2. Then, con- 
figuration Fig. Id is equally probable, and the argument 
can be repeated moving the free spin over large distances 
(spin-charge separation). The final dominant configura- 
tion near a vacancy has the form shown in Fig.lc. The 
key detail to understand the AF enhancement is that in 
the absence of the vacancy the spin at 1 spends roughly 
half the time coupled into a singlet with the spin at 2 and 
the other half with the spin at 0, while with a vacancy 
the spin at 1 forms a singlet most of the time with the 
spin at 2. Then, the 1-2 spin correlation is now enhanced. 
The vacancy has pruned the possible singlet configura- 
tions and the spins next to it no longer "resonate" , but 
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their singlet partners are fixed by geometry. 
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FIG. 1. (a)-(g) Examples of spin singlets relevant for our 
discussion (see text); (h) Cs(l,j) for a S = 1/2 Heisenberg 
model on a 128 sites chain calculated with DMRG. The open 
(full) squares are results with (without) a vacancy at site 0. 
The inset contains the RVB variational results (see text) for 
8-sites with OBC compared against exact results for the same 
chain with OBC and PBC. 1 is the first site of the chain. 

The same pruning of RVB states occurs when correc- 
tions to Fig.la-b are included. Consider a configuration 
with a singlet linking sites 1 and 4 (Fig. If). Such con- 
figuration produces the largest contribution to the spin 
correlation at distance 3. As in the previous case, the 
spin at 1 divides its tendency to form singlets equally 
between sites 4 and -2. But when a vacancy is in- 
troduced, only one possibility for the coupling of the 
spin at 1 with partners at distance 3 remains (Fig.lg), 
enhancing the correlation at that distance. Again the 
nonmagnetic impurity has pruned the singlet configura- 
tions increasing the AF correlations in its vicinity. The 
same argument applies for longer singlets, although mov- 
ing away from the vacancy the AF enhancement is re- 
duced since eventually both Fig.la-b must dominate in 
the bulk. Results for the ground state staggered spin 
correlation Cs(i,j) — (Si • Sj)(— in the standard 
notation, calculated with density matrix renormalization 
group (DMRG) techniques jl5) are shown in Fig.lh for 
the spin 1/2 Heisenberg chain. The AF enhancement is 
large (48% for the bond next to the vacancy). The effect 
can be rephrased as an increase in S(ir) as previously ob- 
served by Eggert and Affleck [[fl] . We here confirm their 
results and provide a simple explanation for the origin of 
the enhancement. 

To provide further support to these ideas, a RVB vari- 
ational calculation on a chain of 8 sites with open bound- 
ary conditions (OBC) was performed. The state used is a 
linear combination of (a) the state with four NN singlets, 
(b) all states with three NN singlets and one singlet of 
length 3, and (c) all states with two triplets of length 2 
and two NN singlets. The relative weights are fixed min- 



imizing the energy. The spin correlations in this varia- 
tional state are also shown in Fig.lh (inset). Contrast- 
ing the results against exact calculations with OBC, it is 
clear that such a simple state contains the main features 
observed in the actual ground state of longer chains. 

The discussion can be easily extended to NN S =1/2 
Heisenberg models with a larger coordination number. 
Consider first a 2-leg ladder. Snapshots of the spin ar- 
rangement in the vicinity of an arbitrary site, e.g. 1, 
would show that this spin spends most of its time singlet 
coupled with spins 2,3 and 4 (Fig. 2a) |ll|. When a va- 
cancy is introduced at site 2, now the spin at 1 forms sin- 
glets only with two partners rather than three (Fig. 2b). 
This pruning of the RVB basis enhances the spin correla- 
tion at distance 1 along the chain, as in ID. Fig. 2d con- 
tains Cs(i,j) for a 2-leg ladder. Shown are the DMRG 
spin correlations both along the leg where the vacancy is, 
and in the opposite leg. The results are in good agree- 
ment with previous calculations |^. The enhancement 
is substantial only for same-leg correlations, and even in 
this case it is smaller than for the S — 1/2 chain, in 
agreement with the picture discussed before (i.e. count- 
ing links the naive enhancement ratio would be 3/2 for 
ladders vs 2/1 for chains). At distance 1 and for same-leg 
correlations, the actual enhancement is ~ 22%. 
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FIG. 2. (a)-(c) Spin singlets relevant for ladders and hy- 
percubic systems (see text); (d) Cs(i,j) for a 2-leg ladder 
calculated with DMRG on a 2 x 32 cluster. The open squares 
(triangles) are spin correlations along the same (opposite) leg 
where the vacancy is located, with a starting site i next to 
the vacancy. The full squares are results without vacancies. 
The vacancy is located at the center of the cluster. 

Numerically the short distance spin correlations 
smoothly interpolate between the results of chains and 
2-leg ladders as the rung coupling J± is reduced with re- 
spect to the chain coupling J. As example, in Fig. 3a 
exact diagonalization (ED) |l~7| results are shown for 
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J±/J = 0.5 using periodic boundary conditions (PBC) 
and OBC. They suggest that the short distance AF en- 
hancement of chains and ladders have a common origin 
independent of the long distance behavior of their spin 
correlations (power-law vs exponential decay). Only the 
range of the Zn-induced disturbance is affected by the 
presence of a spin gap. For example, due to spin-charge 
separation in ID, the S = 1/2 state obtained by removing 
one spin from an even site chain is uniformly distributed 
over the chain ||. However, near the vacancy the AF cor- 
relations are nevertheless enhanced (Fig.lh). Thus, the 
short distance effect discussed here does not seem corre- 
lated with the presence of a S = 1/2 localized state near 
the vacancy. To further confirm these ideas, correlations 
for a J i — J 2 chain ( J 2 being a next to NN coupling) were 
calculated varying J2 / J\ from (no gap in the ground 
state) to the exactly solvable point 0.5 (dimerized gapped 
state). A smooth interpolation between the two regimes 
was observed at short distances (see e.g. Fig. 3b that 
includes a RVB variational calculation). As the Zn im- 
purity density grows, the short distance effects will even- 
tually govern the behavior of the system, and the tail in 
the spin disturbance becomes irrelevant. 

The analysis of models for CuGeC>3 require a special 
discussion. Here a dimerization occurs at low tempera- 
tures H which will be represented by a modulation of the 
orig inal NN exchange using J(l±5)Sj • Sj, with 6 ~ 0.03 
ijbSf . Note that a more proper description of CuGeC>3 
requires the inclusion of phonons to generate the dimer- 
ization dynamically. However, the model used here is 
sufficient for our qualitative studies. The +(— ) sign cor- 
responds to "strong" ("weak") links. The addition of a 
next NN interaction J2/J1 ~ 0.2 — 0.3 is also needed in 
CuGeC>3 |l8|| . Based on the RVB picture discussed be- 
fore it would not be surprising that AF enhancements of 
short distance spin correlations will also appear in Zn- 
doped spin-Peierls chains. In fact we found the surpris- 
ing effect that the enhancement in this case is the largest 
among the family of models studied here. As example, 
consider the results of Fig. 3c, obtained exactly on a fi- 
nite chain with OBC. The pattern of strong vs weak links 
considered here is such that sites 1 and 2 next to the 
end are connected by a weak link. Fig. 3c shows a large 
~ 84% enhancement of Cs(l,2) compared with results 
for a ring. The reason is that for a ring spins 1 and 2 are 
not much correlated since each can form a strong bond 
with a neighboring spin, while for an open chain the spin 
at 1 is "free" and strongly tries to form a singlet with its 
neighbor regardless of the smaller value of the coupling 
,/(l — S) in this link. Thus, the relative enhancement is 
large. If the pattern of weak vs strong bonds is shifted in 
one lattice spacing, then the AF enhancement is smaller 
since the spins of the now strong bond 1-2 are strongly 
correlated for both a ring and a chain with OBC. 
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FIG. 3. C s (i,j) obtained with ED for (a) a 2-leg ladder 
with Jl/J = 0.5 using PBC (full squares) and OBC (open 
squares). The correlations are measured from the edge of a 
2 x 10 cluster; (b) a S = 1/2 chain with J2/J1 = 0.3 and 
14 sites. Open squares denote results from the end of an 
open chain, while full squares are results for a ring. The inset 
shows RVB variational results with the same state used in 
Fig.lh. Open (full) squares are exact results with open (pe- 
riodic) boundary conditions. The triangles denote the varia- 
tional predictions; (c) a dimerized 5 = 1/2 chain with OBC, 
14 sites, and J2/J1 = 0.2, S — 0.03. i is the first spin on 
the chain. The inset contains RVB results on an 8 site open 
chain. The symbols are as in Fig. 3b. The first link in the 
chain near the edge is "weak". 

The local RVB-picture presented here predicts that AF 
enhancement should also occur in dimensions higher than 
1, but with a strength reduced from the results in ID. 
In general, a given spin of an hypercubic lattice would 
tend to form singlets with 2 x D partners in the ab- 
sence of a vacancy and with 2 x D — 1 partners when 
a vacancy is introduced next to it (as exemplified in 
Fig. 2c for D=2). The ideal enhancement factor F would 
be F ~ (2 x D)/(2 x D — 1), but it is certainly reduced 
by the inclusion of corrections beyond NN singlets. In 
Fig. 4a, Cs(i,j) for a (tilted) square cluster of 26 sites 
is shown with and without a Zn-impurity. The relative 
enhancement at distance 1 along an axis away from the 
vacancy is ~ 11%, i.e. smaller than in ID, as predicted. 
Our picture does not depend on the presence of long- 
range order, as exemplified by calculations of Cs{i,j) for 
the frustrated J\ — J2 2D Heisenberg model with a va- 
cancy. Fig. 4b actually shows an AF enhancement near 
the vacancy larger at J2/J1 = 0.5 (where the long-range 
AF order is suppressed) than at J2/J1 = 0.0 (~ 31% 
vs. ~ 11% for the correlation at distance 1). A smooth 
connection between the two results was observed numer- 
ically, and thus once again the local AF enhancement 
seems caused by the RVB character of the short distance 
fluctuations independently of the long distance properties 
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of the model. 




FIG. 4. (a) Cs(i,j) on a 26-sites tilted square cluster calcu- 
lated with ED. Open (full) squares denote results with (with- 
out) a vacancy, i is the site next to the vacancy if present. 
Correlations at distances 1, v2, 2 and \/5 away from the im- 
purity are shown; (b) Same as (a) but introducing frustration 
J2/J1 = 0.5; (c) Alternating contribution to the local spin 
correlations Si (see text and Ref. jnj) for a Ji — J2 — S model 
with J2/J1 = 0.2 and 5 = 0.03 using DMRG on a chain of 160 
sites with OBC (open squares). The full squares are results 
for J2 — 5 = 0.0 on the same chain. In both cases 16 states 
are kept in the DMRG method, and the end of the chain is 
at site 0. 

It is also possible to predict the strength of the en- 
hancement for models with spins S > 1/2. Consider first 
a ID chain. In the limit of S — > 00, the zero temper- 
ature ground state is antiferromagnetically ordered and 
the picture based on RVB singlets is no longer useful. 
In this limit the presence of a nonmagnetic impurity is 
irrelevant for the spin correlations which remain of max- 
imum strength. Assuming a monotonous crossover from 
5 = 1/2 to 5 = 00 a Zn-induced enhancement of Cs(i,j) 
for all S would be expected, but of decreasing magnitude 
as S grows. Indeed for a S — 1 chain the spin correla- 
tion at distance one near an end is ~ 16% larger than in 
the bulk p0| , an enhancement smaller than for S = 1/2 
||19fl . Thus, the vacancy-induced antiferromagnetic en- 
hancement seems to occur in any dimension D and spin 
S, but the effect is the largest for ID S = 1/2 systems 
where quantum fluctuations are crucial. 

The results obtained in this paper allow us to make pre- 
dictions for the NMR spectra of a variety of compounds 
similarly as done before for ID S = 1/2 chains [ pTpfj| ]. In 
Fig. 4c the alternating part of the local spin correlations 
Si = (proportional to the local susceptibility 

Xi) is shown. Here () represents the expectation value 
in the lowest energy state of the subspace of spin one, 



which would contribute to Xi as the temperature or an 
external magnetic field is increased from zero (Si in the 
singlet ground state of a finite lattice would trivially van- 
ish) . The figure shows that both dimerized and undimer- 
ized chains have a nonzero local staggered magnetization 
near vacancies. In the former the enhancement is con- 
centrated in the vicinity of the impurity, while in the 
latter it is spread over all the chain due to spin-charge 
separation in excellent agreement with Ref. pl| . NMR 
experiments should detect a broadening in their spectra 
associated with this effect in several compounds such as 
CuGeOa. At a more speculative level, the AF enhance- 
ment discussed in this paper may be responsible for the 
transition to a 3D Neel order state in Zn-doped dimerized 
chains and ladders Jt],|| ■ Note, however, that this effect 
may depend on the strength of the (small) coupling be- 
tween the structures discussed here, and thus may vary 
appreciably from compound to compound. 
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